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The purpose of BRAHMS is to survey the dynamics of relativistic heavy ion (as weh 
as pp and d-A) colUsions over a very wide range of rapidity and transverse momentum 
[1]. The sum of these data may give us a ghmpse of the initial state of the system, its 
transverse and longitudinal evolution and how the nature of the system changes with time 
[2]. Here I will concentrate on the origin and dynamics of the light flavors, i.e. the creation 
and transport of the up, down and strange quarks. The results presented here are certainly 
not the end of the story. It is my hope that in a few years new detectors will reveal the 
rapidity dependence of the charm and bottom quarks [3]. 

1. Overview of particle production 

Figure [1] shows our preliminary yields of vr^ , and antiprotons versus rapidity for 
central Au+Au collisions at ^/snn = 62.4GeV. The yields of all produced particles are 
well described by Gaussian fits up to rapidities of 0.8 x ybeam- The widths of the two pion 
curves are very close to the prediction by Carruthers for massless particles undergoing 
Landau flow, i.e. cTtt = Inj where 7 is the Lorentz factor of the incoming beams [H El |6]. 
This model assumes that all of the entropy of the system is created at the instant the two 
beams collide and that the system then expands adiabatically until freezeout. One might 
expect that massive particles would have somewhat narrower widths than this, since for a 
given amount of energy transferred they will have lower speeds than lighter particles. The 
data show ordering of the rapidity widths that is the same as at y^s/vw = 200GeV, namely 
(Tfc+ > cj^± > cTfc- > apbar [7J. This is suggestive of a correlation between the production 
of protons and Kaons. We have investigated this link in terms of a chemical analysis of 
quark rapidity densities, as discussed below. 

Since the initial protons are explicity excluded from the Landau/Carruthers model we 
should not expect them to be controlled by the same dynamics. Figure [2] shows our 
distribution of net protons, dN^ /dy — dN^/dy at various energies [8j. These data suggest 
that the average rapidity loss (Ay) of the incoming nucleons rises steadily with the beam 
rapidity up to ^/sNN = 62 GeV before saturating at (Ay) = 2 . 

FigureOshows average values of (mr) versus rapidity for TT^,k^, protons and antiprotons 
for central Au+Au collisions at y^J/vTV = 62.4GeV. The (mr) values are about 10% lower 
than at 200GeV for all particles at all rapidities but the relative drop with rapidity is very 
similar ( 14%). Radial flow gives a uniform velocity boost to all particles and so a bigger 
boost in (my) to the heavier particles. Our data are consistent with a pattern of radial 
flow that weakens as the rapidity increases and/or the beam energy drops. 
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dN/dy vs rapidity in Au+Au 0-5% 
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Figure 1. Preliminary dN/dy versus rapidity from central Au+Au colli- 
sions at ^snn = 62.4GeV 



The fact that the k+ dN/dy distribut ion is wider than the pion distribution may be 
because strangeness is related to the net baryon density. We have investigated this in the 
context of thermal fits to the data. Using the "Thermus package" |14) we have fit five 
independent particle ratios at each rapidity to a grand canonical distribution described by 
a temperature, and two chemical potentials, fig for strange quarks and fig for light quarks. 
More details can be found in [15] . Figure H] shows fig for versus at several different 
rapidities for ^/sJjn = 62.4 and 200 GeV. Note that jj.s ~ l/4//q. The chemical freeze-out 
temperature drops as the rapidity increases or ^snn decreases. For example T = 160 ± 5 
MeV at ^/s]^ = 200 GeV and y=0 and falls to 116 ± 9 MeV at ^/s]^ = 62.4 GeV and 
y=3. 

For non-central heavy ion collisions the geometric eccentricity of the initial state may 
create pressure gradients if the system is strongly interacting. Elliptic flow, i.e. V2 the 2"*^ 
harmonic of the momentum distribution in azimuthal angle, is sensitive to these gradients. 
Since these gradients are self quenching V2 tells us about the early state of the system. 
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Figure 2. Net proton distributions versus rapidity and -^/s/vTV from 
E917,E802 and E877 5GeV [IHl E], NA49 (PbPb) 17GeV [12] and 
BRAHMS 62 and 200GeV ^3\. The 62GeV data are preliminary with 
statistical errors only. 
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Figure 3. Preliminary values of (mr) versus rapidity for tt ,k ,p and p 
for central Au+Au collisions at y^s/vTV = 62.4GeV. The errors are statistical 
only. 

At ^snn = 200 GeV we have studied this phenomena as a function of rapidity for pions, 
Kaons and protons. Figure [5] shows V2 for protons and pions (for clarity the Kaons are 

not shown) as a function of px, and KEt = \Jp't — m as well as V2/nquark versus 

KET/uquark at Central (top) and forward rapidity (bottom). A striking feature of these 
data is how little V2{pt) changes with rapidity. This is despite the fact that when one 
averages over px, (^2) drops by about a factor of two from y=0 to y=3 [16iil7j. The values 
of V2{pt) at all rapidities are close to the maximum values allowed by ideal hydrodynamics, 
under the assumption that all of the initial spatial anisotropy is converted into momentum 
anisotropy. This effect implies very early thermalization. 

Objects produced at forward rapidities will tend to come from the collisions of partons 
with very unequal momenta. Such collisions have a y/s that is lower than the balanced 
collisions which produce particles near mid-rapidity, and so are less likely to produce parti- 
cles at high pt- We find that the mean px of pions, Kaons and protons drops steadily with 
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Figure 4. Preliminary fig versus fig at several different rapidities for 
y^J/Viv = 62.4 and 200 GeV. The errors are statistical only. 



rapidity [7|. Thus even though V2{pt) is rather constant, (^2) drops quickly with rapidity. 
Of course we should not discount the effect of radial and longitudinal flow but this will 
have to wait for a more complete analysis. 

At y=0 the proton and pion data are closest when plotted as V2/nquark versus kT/uguark- 
This effect has been seen by several experiments and suggests that the quarks themselves 
flow until they coalescence into hadrons. This scaling does not work quite so well at y=3 
where V2/nguark at a given kx/ngy^ark is slightly lower for pions than for protons. Since 
only ~ 4% of our reconstructed pions come from k^^^^ decays it is unlikely that this effect 
is caused by a "dilution" of (^2) for the pions. Rather it may reflect a breakdown of 
coalescence, perhaps caused by recombination, in the forward direction. 

The success of non viscous hydrodynamics in describing both the longitudinal and elliptic 
flow suggests that thermalization is archived very rapidly in heavy ion collisions. Black 
holes are a possible generator of very rapid thermalization [18]. Because of the uncertainty 
principle black holes radiate with a temperature that only depends upon their mass, charge 
and angular momentum. For uncharged black holes with zero angular momentum ,T = 
l/(87rGM) where M is the mass, and G is the gravitational constant. The radiation carries 
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Figure 5. Preliminary V2 for protons and tt^ versus pT (left) and KEj- 
(center) and V2/nquark versus KET/uquark (right) for y=0 (top) and y=3 
(bottom) . The data are from semi-central Au-|- Au collisions at y/s^N = 200 
GeV. 



away energy which causes the temperature to increase, further increasing the rate at which 
mass is lost. This effect means that very small black holes would evaporate very quickly, 
converting there energy to a thermal distribution of particles. For charged black holes the 
temperature drops, eventually reaching zero when = GM. One can make an analogue 
of strongly coupled QCD which can be mapped onto general relativity using the ADS/CFT 
correspondence. In this transformation 

M ^ E, i.e. dEx/dy 
Q ^ B, baryon number 
G ^ 1 /2a, string tension 
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and the temperature Tq of a black hole becomes, 



(1) Tq{B)=Tq{Q)^ ^ ' 



(1 + ^1 - 2(tB^/E'^ Y 

which falls very slowly as the ratio B /E increases until 2aB'^ / E"^ reaches 0.9 when it sud- 
denly dives to zero. The temperature Tq(B) is then identified as the chemical temperature 
deduced from an analysis of particle ratios. We make the following working definitions; 

(2) ^= E {^T)dNydy 

i=n,k,p 

(3) B = dNP/dy-dNP/dy 

and take a = 200MeV'^. 

Figures [1] and [3] show that both {nit) and dN/dy decrease with rapidity. At ^/sJiN =200 
GeV both quantities are larger [7J. Therefore dEnp/dy decreases with rapidity and increases 
with ^Jsnn- Figure [2] shows that B = dN^ / dy—dN^ / dy has the opposite dependence, rising 
with rapidity and falling with y/SMN- Thus our largest lever arm for testing Equation [1] 
come from comparing T at y=0 and -^/siVTV = 200GeV with T at y=3 and y^s^viv = 62.4GeV. 
We find that 2aB'^/E'^ changes from (1.0 ± 0.3) * 10"^ to (3.6 ± 0.2) * lO'^ over this range 
while the freeze-out temperature drops from 160ib5 MeV to 116ib9. To accommodate such 
a temperature change within the model would require rescaling our B/E measurements 
by a factor of 5. Using HIJING we estimate that the net baryon number is a factor of 2.0 
times the net proton yield at y=0. At y=3 this factor rises to 2.1. Modifying our definition 
working of B / E to use these estimates and taking vr" = (vr"*" + 7r~)/2 increases B/E by a 
factor of 1.4 ±0.1. 



Conclusions 

The rapidity distributions of particles produced from central Au+Au collisions are 
Gaussian at RHIC energies with widths that increase with ^/sNN. The pion widths are 
well described by the Landau/Caruthers model of a relativistic fiuid expanding at constant 
entropy in one dimension until it hadronizes into massless particles. Of course the motion 
of the particles is not only along the z axis. The increases (m-T) with mass suggests there 
is a radial component to the expansion. The strength of the radial flow decreases with 
energy and rapidity. The fact that the scaled elliptic flow V2{kT) /uquark is proportional 
to kT/uquark at y=0 suggcsts that it is the partons that flow before they coalesce into 
hadrons. There are hints that this proportionality breaks down at y=3 but we need to 
finish our study of the systematic errors of this analysis. The ordering of the rapidity 
widths, (Tfc+ > (T^± > (Tfc- > Gphar docs uot depend upon energy and has been investigated 
in terms of a thermal model of flavor production. The chemical freeze-out temperature 
drops with rapidity and energy while ~ l/4/^<j. All of these data are suggestive of very 
early thermalization. This can be explained by assuming that particles are radiated from 
black holes. We have made the flrst moves towards testing this hypothesis. Currently it 
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seems that this theory predicts too slow a change of the chemical freeze-out temperature 
as the ratio of net baryon number to transverse energy increases. 
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